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Abstract. We study high-density traffic of information packets on sparse modular networks with scale-free 
subgraphs. With different statistical measures we distinguish between the free flow and congested regime 
and point out the role of modules in the jamming transition. We further consider correlations between traffic 
signals collected at each node in the network. The correlation matrix between pairs of signals reflects the 
network modularity in the eigenvalue spectrum and the structure of eigenvectors. The internal structure 
of the modules has an important role in the diffusion dynamics, leading to enhanced correlations between 
the modular hubs, which can not be filtered out by standard methods. Implications for the analysis of real 
networks with unknown modular structure are discussed. 
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1 Introduction 



In recent years network research has been intensified aim- 
I ing at quantitative representation and study of the in- 

■ teractions in complex dynamical systems [1 . These net- 
i works often exhibit hidden structures and inhomogene- 

■ ity at mesoscopic scales. Subgraphs of different sizes and 
topological consistency often appear in real networks, such 
as modules or motifs in gene networks [2j, community 
structure in social networks [3j , topological clusters or dy- 

[ namical aggregation on the Internet 0, and others. It 
has been recognized that in the evolving networks func- 
tional units have emerged, and that in different functional 
networks they may be represented by topologically char- 
acteristic subgraphs, e.g., communities, modules, paths, 
trees, etc. Subgraphs on modular networks can be recog- 
nized topologically by better or tighter connected group 
of nodes [3j . Sparseness of real networks is another feature 
which is tightly connected with the network dynamic sta- 
bility: Large connectivity may induce a chaotic behavior 
(positive Lyapunov exponent) in networks even for sim- 
ple dynamics of its unites [5]. This might be a part of 
the reason why most of the networks in nature (except 
perhaps brain) are sparsely connected [ll|6l[7l[8] . For these 
reasons it is of great importance to understand the inter- 
relationships between dynamics and structure in sparse 
modular networks [9j. This question is also in the focus of 
the present work. We study transport processes on net- 
works with sparsely connected modules by means of the 
numerical simulations of dense traffic, and spectral analy- 



sis of the Laplacian matrix and matrices generated by the 
correlations between the traffic time series. 

Spectral analysis. Important information about com- 
plex network structure and dynamical processes is con- 
tained in the eigenvalue spectra and corresponding eigen- 
vectors of the adjacency matrix and of other, e.g., Lapla- 
cian matrices related to its structure [Tl[T0l[TTl[9]. Recent 
studies of the synchronization of phase-coupled oscillators 
[1^3^12^13] in modular networks have revealed strict rela- 
tionship between the synchronization and smallest nonzero 
eigenvalues of Laplacian matrix [14]. Furthermore, posi- 
tive/negative components of the corresponding eigenvec- 
tors appear to be well localized on the modules [H[T5]. 
Other types of diffusive dynamics, like spreading of dis- 
ease [16] and traffic or navigated walks p!7|[T8|[T9l [2Q] are 
often studied on different networks. The exact relation- 
ship between the return-times distribution of the random 
walk on networks and the spectral properties of the re- 
spective Laplacian matrix, has been derived theoretically 
[TU] and confirmed numerically [9j for the tree graphs. Fur- 
thermore, in [9j the Laplacian spectra have been studied 
in detail for a wide class of sparse and modular networks. 

Time series correlations. In studies of complex net- 
works much effort has been invested in understanding how 
the structure of a network is manifested in network func- 
tion. Beside its theoretical meaning, this question has a 
great practical importance. For instance, in bio-engineering 
[21j and neurosciences [22j, one faces the problem to de- 
sign a network with given function, or to reconstruct a 
network from its measured dynamical output. Often the 
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empirical data are available as time series collected at dif- 
ferent nodes and within specified time windows, for in- 
stance in the market dynamics [23j or gene expressions 
[24] . and other. Different methods have been employed, for 
instance, the correlation matrix reconstruction of neural 
or anti-gene interaction network, [251 126] , gene network re- 
construction algorithms based on SVD and assumed gene 
dynamics model [27^21j. The goal of the reversed engineer- 
ing is to unravel interactions between the nodes which are 
the cause of the observed dynamical output (time series). 
Apart from often limited information, this is a hard prob- 
lem to which both the network structural complexity and 
nonlinear ity of the dynamics contribute. The presence of 
the modules and other mesoscopic inhomogeneities repre- 
senting functional groups of nodes, may increase the com- 
plexity of the problem. In principle, the inherent consis- 
tency of these approaches may rely on certain robustness 
of the inter-dependences between the network structural 
elements and the properties of the time series, seen within 
given type of the dynamics. In particular, time series from 
the diffusion processes, e.g. synchronization in neural net- 
works, are considerably different from the autocatalytic 
regulation in gene networks. The diversity in the activity 
of nodes in the case of random-walk dynamics on struc- 
tured networks is directly related to the node connectivity. 
One of the goals in this paper is to examine the efficiency 
and limitations of the network reconstruction form the 
traffic time series on modular network. For this purpose 
we run known dynamics (traffic of information packets 
with queuing [17]) on known network structure and record 
the traffic time series at all nodes. Then we construct the 
correlation matrix of these time series and use the stan- 
dard filtering methods to uncover the structure behind the 
correlations. The degree of similarity between the filtered 
correlation matrix and the original adjacency matrix, as 
well as between their eigenvalue spectra is quantified. 

The paper is organize as follows: The modular net- 
work structures are introduced and results of simulations 
of traffic of information packets on these networks are pre- 
sented in Section [2l Traffic properties near the jamming 
are studied in detail by statistical means. In Section [3] 
we present the construction and filtering of the correla- 
tion matrix from the traffic time series. Furthermore, a 
detailed spectral analysis of the Laplacian matrices of the 
original networks and the correlation networks is given in 
Section 13. 2[ A short summary and the discussion of the 
results are presented Section [H 

2 Traffic Jamming on Modular Networks 
2.1 Network Structures 

For simulations of the information traffic we use two types 
of modular networks shown in Fig.[Tl in particular: (a) the 
network composed of few large modules with random con- 
nections between them (Net269), and (b) network consist- 
ing of a large number of smaller modules linked through 
a scale-free tree graph. These networks are grown using 
the algorithms which are introduced in [9j . The structural 



properties of these networks are controlled by three pa- 
rameters: the average connectivity M, the probability of 
new module Pq, and the attractivity of node a which con- 
trols rewiring process during the module growth. By the 
numerical implementation and choosing the values of these 
parameters the internal structure of groups (modules) as 
well as the structure of the network connecting different 
modules can be varied in a desired manner (see Ref[9j for 
details). Specifically, for the networks shown in Fig. [T] have 
N = 1000 nodes and the values of the control parameters 
are as follows: M = 2, 10% of rewired links (a = 0.9) 




Fig. 1. (top) Original network with large interconnected scale- 
free modules, Net269, and (bottom) Scale-free tree structure 
with attached scale-free modules. Colors indicate membership 
of nodes to topologically distinct subgraphs. 

and Po = 0.006 for Fig. [D^top), referred as Net269. In 
different structure, shown in Fig. [Tfbottom), we combine 
scale free tree (M = 1, a = 1) with a number of smaller 
modules, where each module is grown taking the same pa- 
rameters M = 2, a = 0.9. The modules are attached to 
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ending nodes of tfie scale free tree Fig. [TJbottoni). Note 
tfiat tfie internal structure of tfiese modules is statistically 
tfie same as in case of Net 269, fiowever tfieir size varies be- 
tween 20 and 50 nodes and tfieir number is cfiosen in sucfi 
a way tfiat approximatefy fiaff of totaf number of nodes are 
members of tfie modufes, wfiife tfie otfier fiaff of nodes be- 
fong to tfie underfying tree structure. It sfioufd be stressed 
tfiat, in contrast to tfie Net269 in Fig.IlJtop), wfiicfi can be 
fuffy partitioned into modufes (communities), tfie modu- 
far network in Fig. [Tfbottom) consists of modufes and tfie 
underfying network (tree), wfiicfi is structured itseff. In 
tfie foffowing sections we wiff sfiow fiow tfiese two compo- 
nents of tfie modufar networks affect tfie traffic and study 
tfie patterns of correfations in tfie traffic signafs. 

2.2 Traffic of Information Packets 

In tfiis Section we simufate dense traffic of information 
packets on tfie network structures introduced above. We 
use our traffic modef introduced earfier (see review articfe 
[17] for detaifs). We describe in sfiort tfie main features 
of tfie traffic modef and tfie refevant parameters (detaifs 
of tfie numericaf impfementation of tfie rufes are given in 

my- 

— Creation and assignment. At eacfi time step eacfi node 
creates a packet witfi a given rate R and assignees it a 
randomfy sefected recipient node (defivery address). 

— Navigation. Eacfi node processes a packet from top 
of its queue (LIFO-queue) towards one of its neigfi- 
bours. Tfie neigfibour node is sefected according to 
nnn-navigation rufe [29t fT7] , in wfiicfi tfie node searcfies 
for tfie packet's recipient address in its neigfibour fiood 
witfiin two-fayer deptfi. If tfie recipient is not in tfie 
searcfied area, tfie packet is sent to a random neigfi- 
bour, wfio repeats tfie searcfi in its neigfibourfiood and 
so on. 

— Queuing. Wfien more tfian one packet is found at tfie 
same node, tfie packets make a queue in tfie buffer 
at tfiat node, waiting to be processed. We use a fixed 
maximum buffer size H = 1000 packets at each node. 
If tfie buffer of a sefected node is fuff, i.e., as at tfie 
jamming thresfiofd, tfie packet can not be defivered 
and waits for a further possibifity to be defivered. One 
packet per time step is processed. 

— Delivery. When the packet arrives to its destination 
(recipient node) it is removed from the network. 

Simufations reveafs that the traffic properties depend on 
the parameters, i.e., posting rate maximum queue fength 
the queuing discipfine (LIFO), and the search depth 
{d = 2 in the case of nnn-search), as weff as on the struc- 
ture of the underfying network. Note that the transport of 
packets with the nnn-search differs from the random diffu- 
sion (random- wafk dynamics) in that, when the recipient 
node is found within next-neighbourhood of a processing 
node, the packet goes directfy to the recipient node. Oth- 
erwise, the packet performs a random wafk. In this way, 
the traffic on networks with structure which is more suit- 
abfe for the nnn-search is much more effective. In particu- 
far, it was shown in Refs. [30j that the cfustered scafe-free 



network of the structure of Webgraph (i.e., with two hub 
nodes and a farge number of triangfes) is most suitabfe for 
next-neighbour-search navigation. The statisticaf features 
of the traffic on the Webgraph structure (with parame- 
ter a = 0.25) compared to the scafe-free tree (see detaifs 
in [29]) show, among other measures, shorter travef and 
waiting times of packets, and up to 40 times farger packet 
density before the jamming point is reached. These fea- 
tures are important for understanding the traffic on our 
modufar networks shown in Fig. [U where, as expfained 
above, the internaf structure of each modufe in Net269 
Fig. mtop) has the Webgraph structure (with a = 0.9) 
and the supporting network in Fig. (ifbottom) is a scafe- 
free tree. 

In Fig. [2ja) bottom panef, we show severaf resufts 
of the traffic on the modufar network Net 269. Measured 
are local and global statisticaf properties for fixed post- 
ing rate R = 0.8. The probabifity distributions of travef 
times of packets P{Tt) between creation and defivery and 
the distribution of waiting times of packets at different 
nodes afong the path, P{tw) are broad range distributions, 
which are characteristic for the networks structure and 
given packet density (befow jamming point). Specificaffy, 
the waiting time distribution shows the power-faw decay 
P{tw) ^ ^ with sfope ^ 2.1, indicating that that 
at this posting rate R the network operates cfose to the 
threshofd of jamming (divergence of the average waiting 
time, which is compatibfe with r^^ < 2, is one of strik- 
ing features of the traffic jamming). The distribution of 
travef times appears to be affected with the network mod- 
ufarity. In contrast to the Levy-type distribution on the 
Webgraph (see [29j), the modufes of the same structure in- 
terconnected as in Fig.lljtop), fead to a distribution with 
weak power-faw at smaff times, and a q-exponentiaf taif 
(see recent work [3T] and references therein): 

P{X) = AX-^ X (1 + (1 - g)X/Xo)^/(^-^) ; (1) 

with r ~ 0.33 and q = 1.29. A simifar expression with two 
sfopes (r ~ 1 and q ~ 1.21) can fit the distribution P{At) 
of time intervafs between successive events at a node (fit 
not shown). Note that in the case of dense traffic that we 
simufate here, the distribution P{At) refers to return of 
the activity (not the packet) to the same node, and thus 
it is different from the famifiar autocorrefator, the return 
time of the random wafk to the origin (see [9j for the simu- 
fations of the autocorrefator). The appearance of two dif- 
ferent sfopes in the distribution of return of the activity to 
the node P{At) in Fig.[2](a), suggests uneven rofe of differ- 
ent nodes in the traffic. This is further studied in terms of 
number of packets {hi{tk)} processed by each node within 
a given time window Twin = 1000 time steps, shown in 
Fig. EJb). The time series {hi{tk)}^ where tk stands for 
the index of successive time windows, is recorded at each 
node of the network i = 1, 2, • • • TV. Dispersion of each of 
these time series, a[i] is pfotted in Fig.[2l^b) bottom panef, 
against its average vafue < h[i] >. In this scatter pfot each 
point stands for one node of the network. The pfot shows 
the scafing behavior 

a\i] = constx < h\i] >^ ; (2) 
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<h[i]> 



(a) (b) 

Fig. 2. (a) Distributions of ttie travel times (Ptt) and waiting times (Ptw) of packets, (bottom panel), and the distribution 
of time intervals between successive activity of a node averaged over the network, (top panel), (b) Scatter plot dispersion a[i] 
vs. average < h[i] > of the time series of all nodes, (bottom), and the ranking distribution of nodes according to the average 
number of packets < h[i] > (o) and the distribution of node occupation within time window Twin = 10^ steps (x). All data on 
the Net269 (Fig. [Hop) for posting rate R = 0.8. 



with /i G [1/2,1], which is characteristic for many real 
dynamical systems (see recent review [32 ). 

The origin of such scaling in traffic models on complex 
networks has been attributed [33j to node groupings ac- 
cording to their role in the traffic, which, in turn, is related 
to their topological (or dynamical) centrality. At the tip of 
the plot in Fig. EJb) bottom, one can distinguish a group 
of 10 most active nodes. The same group is also separated 
in the ranking plot (flat part at the beginning of the curve 
in the upper panel) from the rest of nodes. The ranking 
plot exhibits a power-law decay (Zipf's law) according to 

< h[i] >~ rank~^ ; (3) 

with 7^1, where ranking ranki of nodes i = 1, 2, • • • A/" 
is made according to the average number of packets < 
h[i] > processed by the node i. Furthermore uneven role 
of nodes is demonstrated in the probability distribution 
of the number of packets processed within the time win- 
dow < h[i] > (or the "occupation probability" of a node 
p[i] =< h[i] > /Twin)' The distribution is also shown in 
Fig. [2jb) top panel. The central peak appears resembling 
the "ergodic" system behavior [34], however, at the side 
of large occupancy additional peaks are present represent- 
ing the activity of modular hubs and other highly active 
nodes within the modules. Note that for this traffic den- 
sity (posting rate) there the nodes which have maximum 
occupancy p = 1 may occur, which is another signature of 
the pre-jamming behavior in the network. 

The activity of different nodes and their role in the 
traffic process on the modular network Net269 are shown 
comparatively in Fig. [SJa). The most active nodes (red) 



are modular hubs. Gradually lesser activity is visible at 
nodes within each module, according to their reduced cen- 
trality (scale-freeness of the internal module structure). 
For the comparison, the 3-dimensional plot of the time se- 
ries of all nodes in the scale-free tree with attached mod- 
ules is given in Fig.[3l^b). Jamming on the tree graph (indi- 
cated in red color) is visible at the place of modules. This 
indicates that packets often get confined within a mod- 
ule due to ineffective search mechanism on the underlying 
tree graph. In the following section we analyse the corre- 
lations between these time series and the spectra of the 
correlation matrix. 



3 Correlations and Spectra 

3.1 Correlation of traffic time series 

The elements of the correlation matrix Cij are obtained by 
calculating the Pearson's correlation coefficient between 
time signals hiitk) for each pair z,j of nodes in the net- 
work, given by: 

^ _ E^J^^fa)- < hi >][hj{tk)- < hj >] 



Here hi{tk) is activity of node i in time- window t^^ < hi > 
is average activity during the whole time period, 

r7+ ^ — ^ 
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Time window Time window 

(a) (b) 

Fig. 3. Tliree- dimensional plot of the time series of the number of packets processed by a node within a time window Twin — 10^ 
steps, {/li(tfc)}, for all nodes in the network i = 1, 2, • • • 1000 and for 132 time windows for Net269 (a) and for tree-with-attached 
modules (b). Color map: dark-blue corresponds to low, while dark- red to highest recorded value. 



where rit is the number of time windows considered, and 
(Ti is standard deviation of the time signal at node i: 

= • 6 

nt 

The Pearson correlation coefficient takes values from —1 
(strong anti-correlations) to +1 (strong correlations be- 
tween nodes). Using Eq. (|4]) we obtain the correlations 
between nodes in modular networks from the time signals 
obtained in simulations of traffic, presented in Fig. [3l For 
this kind of time signals, the distribution of the correlation 
coefficients P{Cij) strongly depends on the overall traffic 
density (or posting rate R), and for the posting rate ap- 
proaching the jamming threshold used in this simulations, 
the peak of the distribution is moved towards right edge. 
Specifically, for the signals shown in Fig. [3l the values of 
correlations are centered around c = 0.2 (Net269) as to 
c = 0.35 (tree with attached modules), with negligible 
density of the negative correlations. Nodes connected in 
the original network have high positive correlation coeffi- 
cient which is a result of their similar activities during the 
measurement time. The 3-dimensional plot of the correla- 
tion matrix for the traffic signals on the modular network 
Net269 is shown in Fig. [D^a). Similarly, the correlation 
matrix of the signals recorded on the tree network with 
attached modules is shown in Fig. [SJtop). It is clear that 
in both networks, the most active nodes in each module 
(module hubs) have a high correlation coefficient with each 
other and with the rest of the nodes (cf. Figs. [Hand [5]), 
although topological connections between them might not 
be present. 

In order to extract the information about network struc- 
ture from the correlation matrix, we first observe that only 
correlations above certain threshold Co, i.e., Cij > Co 



might be relevant. In this way one attempts to separate 
the values of potentially relevant correlations from those 
which arise accidentally (and are normally distributed around 
central peak). Although many weak correlations are fil- 
tered out in this way, the remaining matrix still contains 
many spurious links, compared with the original adjacency 
matrix of the Net 269. This is demonstrated graphically in 
Fig. m^a), where the correlation matrix of the traffic sig- 
nals is shown with threshold Co = 0.4. Note that with 
this threshold the remaining correlation matrix contains a 
single connected component. Higher threshold values may 
result in fragmentation. 

As the Fig.lD^a) shows, the correlation matrix contains 
already information about modules (diagonal blocks) and 
their size. The picture is much less clear in the case of small 
modules on the scale- free tree (cf. Fig. [5] (top)). Generally, 
the hub of the scale-free tree has strong correlations with 
other nodes, because of the large signal (large number of 
packets processed) on the hub: first row and first column 
in the matrix. Similarly, the correlations are enhanced be- 
tween the hubs of the large modules in Fig. [H^a), since 
the hubs inside the modules carry the largest traffic. For 
small modules on the tree, the walker gets trapped inside 
the module for longer time, since the module is linked 
to the rest of the network by a single node. Hence, en- 
hanced correlations between the modules remain above 
the threshold and are seen as the off-diagonal blocks in 
Figs, m and m Generally, more sophisticated methods are 
necessary in order to reduce the number of such spurious 
correlations, which are not related with the occurrence 
of a direct link between the nodes in the adjacency ma- 
trix [351 125]. Here we apply one of the filtering methods 
which utilizes the affinity transformation [25 | [26 ] . The idea 
is to multiply each element Cij of the correlation matrix 
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Fig. 4. Correlation matrix of the traffic signals on (a) Net269 and (b) on same network with randomized links inside the modules, 
(c) and (d): Corresponding filtered correlation matrices. In all cases shown are links above the threshold value Co = 0.4. 



with a factor Mij which is constructed from the elements 
the rows i and j in the correlation matrix in the follow- 
ing way: Excluding the diagonal elements Ca and Cjj the 
remaining matrix elements of the correlation matrix are 
first reordered to form the n = (TV — l)-dimensional vec- 
tors {Cij.Cii, ...,Cin} and {Cji, Cji, ...,Cjn}- Then Mij is 
computed as the Pearson's coefficient of the components of 
these vectors. Then the matrix element of the filtered 
correlation matrix is given by the product 

= Mi^Cij . (7) 

In this way, the correlation between the nodes i and j is en- 
hanced if the corresponding meta-correlation element Mij 
is large (i.e., the nodes i and j see the rest of the network 
in a similar way), and reduced otherwise. If two nodes are 
linked on the network, it is expected that their correla- 
tions with other nodes are similar, resulting in value of 
their meta-correlation close to one, otherwise coefficient 
M is closer to zero. Hence, the multiplication of the el- 
ements of the correlation matrix with meta-correlations 
should increase the difference between true and random 
correlations. For the filtering procedure we use the whole 
correlation matrix without any threshold. Note also that, 
in contrast to most general case where a shift of the inter- 
val [-1,1] to [0,1] is performed before the filtering, our cor- 



relation matrix of traffic signals is naturally shifted to the 
positive side. The filtered correlation matrix is also shown 
in Fig. [41(c), only the links stronger than the threshold 
are shown. In comparison with the unfiltered correlation 
matrix, the number of matrix elements is considerably re- 
duced relative to the same threshold value. The effects are 
much more prominent in the case of the tree with small 
modules, shown in Fig. [5] (bottom). 

We further notice that the internal structure of the 
(large) modules play a role in the dynamics and the corre- 
lation patterns. The occurrence of hubs, better connected 
nodes which carry most of the traffic inside the modules, 
leads to strong correlations of the traffic signals with other 
nodes within the module and between different modules. 
Some of these correlations can not be filtered out, as seen 
in Fig.m^c). In order to support this conclusion, we made 
random rewiring of the links inside each module by keep- 
ing the total number of links and links between the mod- 
ules unchanged. Then we run the traffic on the randomized 
network and construct the correlation matrix of the traffic 
signals. The results are shown in Fig. [H (b) and (d): the 
correlation matrix of the randomized network and the cor- 
responding filtered correlation matrix (lower panel), sug- 
gesting that the filtering procedure is much more effective 
in the case of randomized structure of the modules. 



Bosiljka Tadic, Marija Mitrovic: Traffic Jamming & Correlations Patterns 



7 




P 200 iOC EIU IN 1«H 



Fig. 5. Correlation matrix (top) and filtered correlation ma- 
trix (bottom) obtained from the traffic time series on tree with 
attached modules. Shown are the links above the threshold 
value Co = 0.4 . 



By applying the threshold in the filtered correlation 
matrix one can visualize the graph structure (shown in 
Fig. [6l top) for both original and randomized version of 
modules. Compared to the original adjacency matrix, these 
structures contain many spurious links, although the num- 
ber of such links is considerably smaller in the randomized 
version. It is interesting to note that the association of the 
nodes to given modules is almost entirely preserved as in 
the original network. In Fig. [6] the links represent the fil- 
tered correlation matrices, but each node carry the color 
which indicates its membership to a module in the orig- 
inal network Net 269. Again, the randomized connection 
inside the modules lead to the correlation network with 
clear modular structure with only few nodes wrongly at- 
tributed to their original modules. 



3.2 Spectral analysis 

Further details of the interdependences between the diffu- 
sion on networks and their structure are obtained through 
the spectral analysis of the Laplacian operator. The de- 
tailed spectral analysis of the normalized Laplacian re- 
lated to the modular networks has been reported in Ref. 
[9]. The spectral density and the structure of the corre- 
sponding eigenvectors show specific features which are re- 
lated to the modularity and other properties (clustering, 
average connectivity, etc.) of the networks. 

Here we perform spectral analysis of the normalized 
Laplacian related to the correlation matrices C and 
for both types of the modular network structures discussed 
above. In order two exclude selfedges, all elements on di- 
agonal are set up to zero value. The normalized Laplacian 
related to random- walk type dynamics is given by [QIITQ] 

A 

Lij = dij ; , (8) 

where Aij are elements of adjacency matrix of the graph, 
and Qi and qj are the degree of node i and j. It has lim- 
ited spectrum in the range [0, 2] and an orthogonal set of 
eigenvectors, which makes it suitable for numerical study 
and comparisons of different structures. For the networks 
extracted from the correlation matrix C and filtered corre- 
lation matrix the Laplacian is obtained by construct- 
ing the binary graph A in Eq. ([8]) with the elements of the 
matrix C or C^, where Cij > Co are set to unity, and 
zero otherwise. As we show in our previous work [9j, the 
spectral properties of the Laplacian depend on the net- 
work topology, and can be used for the identification of 
its mesoscopic structure. 

The spectrum of the Laplacian related to the origi- 
nal Net269 has six lowest nonzero eigenvalues. Fig. [6](c), 
which are separated from the rest of the spectrum, and 
the largest eigenvalue Xmax < 2. The number of smallest 
eigenvalues (A > 0) are correlated with number of well 
separated subgraphs in the network, while the number of 
eigenvalues A = corresponds to number of disconnected 
components [9]. The spectrum of the networks obtained 
from C, also shown in Fig. [BJc), has 15 zero eigenvalues 
or disconnected components for the applied threshold Co- 
The maximal value A^^ = 2 also indicate that a sub- 
graph with tree-like structure (chain of nodes) [91136] oc- 
curs, in contrast to the original network Net269. The num- 
ber of single nodes or disjointed subgraphs containing a 
few nodes in the network increase after filtering, see Fig. 
m Furthermore, meta-correlation network has continuous 
spectrum up to zero, although the network exhibits mod- 
ular structure. The absence of the gap between the lowest 
eigenvalues and the rest of the spectrum, which is char- 
acteristic for modular networks, is due to many spurious 
links between hubs which increase connectivity between 
modules, as also seen in Fig. IH^a). The high density of 
links, and thus higher average connectivity in correlation 
networks, affects the spectrum in the middle part around 
A^ = 1, shown in Fig. [6jc). In the case of the network 
with randomized links in the modules, shown in Fig.[6l(b), 
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Fig. 6. Weiglited networks obtained from filtered correlation matrices for traffic on Net269 (a) and on network with randomized 
module structure (b). Shown are links above the threshold Co = 0.4. Colors indicate the original membership of nodes to different 
modules of the Net 269. Ranking of eigenvalues of the Laplacian matrix for networks obtained from correlation matrix and filtered 
correlation matrix compared to one obtained for original network Net269 (c) and randomized network (d). 



the situation is more clear: the spectra with the filtered 
correlation matrix and the real adjacency matrix coincide 
at both ends, which is compatible with the precise struc- 
ture of the modules in the top panel of the same Figure. 
In the intermediate part, however, the deviation between 
the filtered and real adjacency matrix are large, suggest- 
ing that this part of the spectrum can not be effectively 
used for the identification of the true network structure. 
In the case of the scale-free tree with attached small mod- 
ules the analysis of the spectra of the filtered correlation 
matrix is entirely ineffective (not shown), suggesting that 
other methods are necessary for the identification of their 
modular structure from the dynamical time series. 

Spacing distribution is another feature the spectra which 
is often studied in the context of the random matrix the- 
ory. Recently it was shown that the distribution of nearest- 
neighbor spacing of the eigenvalues for various model net- 
works (i.e., small world, random networks and scale free 
networks), follow the universal Gaussian orthogonal en- 



semble (GOE) of random matrix theory [37] • In this part 
of the paper we study the spacing of the eigenvalues in 
various modular networks. We show that, although the 
spectra of the modular networks exhibit characteristic fea- 
tures related to their mesoscopic structure, connectivity 
and clustering (see [9j), the spacing distribution of their 
spectra can not be distinguished from the ones in the ran- 
dom matrices. In Fig. [71 the distributions for the spacing 
of the eigenvalues of the adjacency matrices of different 
modular networks with varied structural parameters are 
shown. We calculate nearest neighbor spacing distribution 
for modular networks with different average connectivity, 
M, and different clustering coefficient, grown according to 
the algorithms in [9j, and compare it with non-clustered, 
non- modular scale free networks. 

We use the standard procedure to determine the nearest- 
neighbour eigenvalue spacing (see [37J and references therein). 
The eigenvalues of the adjacency matrix are first ordered 
as Ai < A2 < . . . < Aat, where N is the size of the network. 
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Fig. 7. Distribution of eigenvalues spacing Si = A^+i — 
for non- modular network M = 10, a = 1 and for modular 
networks with different average connectivity and different value 
of clustering coefficient: Net269 (M = 2,a = 0.9, Po = 0.006), 
Net569 (M = 5, a = 0.9, Po 0.006) and Net566 (M = 5, a = 
0.6, Po = 0.006) 



Then the "unfolding" procedure is performed, which can 
be accomplished using the transformations = Nav{^)^ 
where Nav{^) is the cumulative distribution of the eigen- 
values. For a discrete spectrum, the transformation en- 
ables to extract the "smooth" part from the sample spec- 
trum and separate it the "staircase" in the cumulative 
distribution. The analytical form of the distribution Nav 
is obtained using polynomial fit of "smooth" part. The 
spacing distribution P{s) for the random-matrix models 
is then defined as the probability of finding the nearest- 
neighbor eigenvalue in the spectrum at the distance s. 
For the unfolded eigenvalues the distance is defined as 
Si — ^i+i — ^i. In the random matrix theory P(s) is de- 
scribed by Brody distribution 



P = As^e"^''^^ , (9) 



where A = (1 + /?) and a = [r(f^)]^+^ Brody distri- 
bution is an empirical formula characterized by one pa- 
rameter f5. For the case of Gaussian random matrices [38j 
and various model networks [37j it appears that (5 ^ \. 
We perform the same procedure described above for four 
types of networks with six modules and with different clus- 
tering and average connectivity. The spacing distributions 
are shown in Fig|7l compared to the analytic form of the 
Brody distribution in Eq. ([9]) with /3 = 1. The results 
show that within the numerical error, even though clus- 
tering, modularity and average connectivity influence the 
spectral density of the networks, the spacing distribution 
remains unchanged and can be fitted with the same pa- 
rameter /3 ^ 1. 



4 Conclusion 

We have simulated high density traffic of information pack- 
ets with local search and queuing at nodes [17J on two 
types of networks with higher structures — modules of dif- 
ferent sizes and scale-free internal structure. One of the 
goals was to emphasize the role of modular structure against 
the structure of the underlying network connecting these 
modules. Our results suggest that the network composed 
of several modules with clustered scale-free structure can 
bear much larger traffic density before jamming occurs, 
compared to the network of the same size but with a sin- 
gle module of the same structure. However, the traffic ef- 
ficiency, measured with the statistical parameters and the 
scaling exponents, is reduced and strongly dependent on 
the way the modules are interconnected. Particularly, the 
slope of the travel time distribution is reflecting the struc- 
ture of the connective network. Whereas, the traffic jam- 
ming first occurs in the modules due to the trapping of the 
packets within a module away from their destinations. 

Another aspect of this work concerns the reversed prob- 
lem: recovering the network structure from the correla- 
tions between the traffic time series. In our approach we 
generate the network-wide time series (traffic signals at 
each node) on pre-defined network structures, and have 
demonstrated how some standard filtering procedures works 
in the presence of modules. Particularly, the internal struc- 
ture of the modules (presence of hubs and hierarchy be- 
tween the nodes) induces spurious correlations which are 
elusive for the filtering methods. Our results indicated 
that additional input about how the modules might be 
structured is necessary in order to increase the validity of 
the filtering methods. Two cases should be differentiated: 
One, when the network is fully partitioned into subgraphs 
(modules) of similar structure, like our Net269, and the 
other, when the modules are immersed into connecting 
network which has a nontrivial structure by itself. Our 
findings are corroborated with the analysis of the eigen- 
value spectra of Laplacian for both the original and the 
correlation (filtered) matrices. Specifically, we have con- 
sidered the spectra of the normalized Laplacian of these 
matrices, which has the bounded spectrum in the range 
[0,2] and it is suitable for the comparison: When the fil- 
tering is improved, the spectra of the filtered correlation 
matrix converge towards the spectra of the real connec- 
tivity matrix. The convergence first occurs in the edges of 
the spectrum, and the eigenvectors related with the low- 
est nonzero eigenvalues localize on the modules, although 
the number of links between these modules can still be 
unrealistically large. We have shown that the efficiency 
of the procedure is strongly dependent on internal homo- 
geneity of these modules. More homogeneous modules can 
be detected with higher accuracy. These results may have 
implications for the real complex systems, such as gene ex- 
pressions [21,24^26j, or stock market data [39], where the 
time series are used to reconstruct the (unknown) under- 
lying network structure. Our message is that, in contrast 
to local node connectivity, the mesoscopic structure (func- 
tional modules) can be identified with better accuracy. 
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